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We study the impact of a magnetic field, generated in collisions of relativistic heavy ions, 
on the decay probability of a quarkonium produced in the central rapidity region. The quark 
and antiquark components are subject to mutually orthogonal electric and magnetic fields 
in the quarkonium comoving frame. In the presence of an electric field quarkonium has a 
finite dissociation probability. We use the WKB approximation to derive the dissociation 
probability. We found that quarkonium dissociation energy, i.e. the binding energy at which 
the dissociation probability is of order unity, increases with the magnetic field strength. It 
also increases with quarkonium momentum in the laboratory frame due to Lorentz boost of 
electric field in the comoving frame. We argue that J/ip's produced in heavy-ion collisions at 
LHC with P± > 9 GeV would dissociate even in vacuum. In plasma, J/ip dissociation in a 
magnetic field is much stronger due to the decrease of its binding energy with temperature. 
We discuss phenomenological implications of our results. 



I. INTRODUCTION 



In their seminal 1986 paper [1] Matsui and Satz suggested that quarkonium production in 
high energy nucleus-nucleus collisions can be used as an indicator for production of quark-gluon 
plasma (QGP). They argued that a quarkonium of radius r will dissociate in QGP when the Debye 
screening radius rn{T) becomes smaller then r. The observed effect would be an "anomalous" 
suppression of quarkonium yield. In practice, determination of the screening radius r£)(T) turned 
out to be a remarkably difficult problem even in a static medium. The existing approaches to solve 
this problem include lattice QCD calculations of quarkonium correlators [2-7], construction of 
potential models of quarkonium spectral functions [8-14], and use of effective field theory [15-17]. 
It is remarkable that in spite of much progress there still exists substantial uncertainty in the value 
of the J/tp dissociation temperature and in the functional form of ro(T), see e.g. [18, 19]. Another 
complication arises due to "cold nuclear matter effects", i.e. nuclear effects that are independent 
of temperature and that would occur even if no plasma were formed. There is ongoing controversy 
as to the nature of these effects, although there is agreement that they lead to suppression of 
quarkonium yield [20-29]. 

It has been recently realized [30] that colliding heavy-ions produce a very strong magnetic field 
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in the direction perpendicular to the reaction plane (defined as a plane containing momenta of the 
ions and the impact parameter) . This has a number of interesting phenomenological consequences 
[30, 32-34]. The strength of the magnetic field at RHIC is estimated as cBq ps m^, and at LHC 
as cBq « 15 777,2 [30, 31]. The relaxation time of the magnetic field neglecting the back reaction 
is of the order of the width of the two Lorentz contracted nuclei divided by c, which is ~ 0.1 fm 
at RHIC and ~ 0.01 fm at LHC. However, we argued recently [32] that the relaxation time of 
this field is actually much larger if the back reaction is taken into account.* In fact, the magnetic 
field can be treated as static if the distance over which it significantly varies is much larger than 
the quarkonium radius. We will refer to this approximation as "quasi-static". For a quarkonium 
with binding energy and radius a s /eb, the quasi-static approximation applies when e^t ja s S> 1. 
For r = 2 fm we get e^t ja s ps 23, which is comfortably large enough to justify the quasi-static 
approximation, where we assumed that e& is given by its vacuum value. As T increases e& drops, 
while r increases. Temperature dependence of e& is model dependent, however it is certain that 
eventually it vanishes at some finite temperature To. Therefore, only in the close vicinity of To, 
i.e. at very small binding energies, the quasi-static approximation is not applicable. We thus rely 
on the quasi-static approximation throughout the paper. 

A magnetic field has a three-fold effect on quarkonium: 

1. Lorentz ionization. Suppose quarkonium travels with constant velocity in a magnetic field 
in the laboratory frame. Boosting to the quarkonium comoving frame, we find mutually 
orthogonal electric and magnetic fields given by Eqs. (3), (4). In the presence of an electric 
field quark and antiquark have a finite probability to tunnel through the potential barrier 

* This can be seen as follows. Upon collision, the magnetic field Bo rapidly decreases with time inducing an electric 
field of strength E ~ BqR/t circulating around the direction of B ; here R is the nuclear radius and r the 
relaxation time. This electric field generates circular currents that, by Lenz law, support the original magnetic 
field. There are two types of generated currents: (i) Foucault currents in the QGP; (ii) Current of charged fermions 
produced via the Schwinger mechanism [35]. Note that lepton Schwinger pairs are certainly not in equilibrium 
with the QGP and therefore do not contribute to the plasma electrical conductivity. In the first case, for a medium 
with electrical conductivity a we estimate Bo ~ jR ~ aER ~ aBoR 2 /r, which implies 

r ~ aR 2 . (1) 

Using the values of the electrical conductivities obtained in [36] and in [37, 38] we get r ~ 8 fm and r ~ 0.4 fm 
correspondingly. In the second case, recalling that the density of the Schwinger pairs is n ~ (eE) 2 T we have 
Bo ~ envR ~ e 3 B 2 R 3 /r, which implies 

r~e 2 (eB )B 3 . (2) 

We took into account that the Schwinger pairs are relativistic v ~ 1 since their longitudinal momentum is of order 
eEr ~ eBoR ^> m e for R — 5fm. We estimate that at RHIC r ~ 0.5 fm, while at LHC r ~ 8 fm. A more accurate 
estimate of r requires not only better knowledge of electrical conductivity in expanding plasma, but also numerical 
solution of the magneto-hydrodynamic equations. 
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thereby causing quarkonium dissociation. In atomic physics such a process is referred to 
as Lorentz ionization. In the non-relativistic approximation, the tunneling probability is of 
order unity when the electric field E in the comoving frame satisfies eE > m 1 ! 2 ^ 2 (for 
weakly bound states), where sj, is binding energy and m is quark mass, see (26). This effect 
causes a significant increase in quarkonium dissociation energy Ed- The corresponding results 
for J/ip are exhibited in Fig. 4. They suggest that J/tjj in plasma dissociates at a much lower 
temperature then it would in the absence of a magnetic field. 

2. Zeeman effect. Energy of a quarkonium state depends on spin S, orbital angular momentum 
L, and total angular momentum J. In a magnetic field these states split; the splitting energy 
in a weak field is AM = ^gJ z , where J z = — J, — J + 1, . . . , J is projection of the total 
angular momentum on the direction of magnetic field, m is quark mass and g is Lande factor 
depending on J, L and S in a well-known way, see e.g. [39]. For example, J/ip with S = 1, 
L = and J = 1 (g « 2) splits into three states with J z = ±1,0 and with mass difference 
AM = 0.15 GeV, where we used eB$ = 15m 2 (the value of magnetic field expected at 
LHC [30, 31]). Thus, the Zeeman effect leads to the emergence of new quarkonium states in 
plasma. 

3. Distortion of the quarkonium potential in magnetic field. This effect arises in higher order 
perturbation theory and becomes important at field strengths of order B ~ 37rm 2 /e 3 [40]. 
This is 3n/a times stronger than the critical Schwinger's field. Therefore, this effect can be 
neglected at the present collider energies. 

The purpose of this paper is the analysis of the Lorentz ionization of quarkonia. 

Our paper is organized as follows. In Sec. II we give fully relativistic derivation of the Lorentz 
ionization probability w using the WKB approximation. We treat the force binding q and q into 
quarkonium as a short-range one i.e. (Me^^R <C 1, where e& and M are binding energy and 
mass of quarkonium, respectively, and R is the radius of the nuclear force given by R ps (as/a) 1 / 2 , 
where a = 1 GeV/fm is the string tension. For example, the binding energy of c and c in J/ijj 
in vacuum is = 0.64 GeV <C M/R 2 = Ma/a s ~ 3 GeV. This approximation is even better 
at finite temperature on account of e& decrease. Regarding J/ip as being bound by a short-range 
force enables us to calculate the dissociation probability w with exponential accuracy w ps e~f , 
independently of the precise form of the quarkonium wave function. This is especially important 
since solutions of the relativistic two-body problem for quarkonium are not readily available. We 
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will argue that exponential accuracy is adequate for determining the quarkonium ionization energy 
due to rapid variation of / with binding energy. Our result for w is in agreement with previous 
calculations [41]. 

In Sec. Ill we study non-relativistic approximation of ionization probability w. We argue that 
it provides a remarkably good estimate of relativistic formulas, see Fig. 2. This is an important 
observation as it allows us to include the contribution of the quark spin interaction with the 
magnetic field. A fully relativistic calculation that accounts for the spin contribution is not yet 
available. This is discussed in Sec. HID. In Sec. Ill we also show that relativistic formulas for w 
derived in Sec. II reduce to non-relativistic formulas found in the literature. 

In Sec. IV we calculate the dissociation energy of J/ip m. a, magnetic field and claim that it 
strongly depends on magnetic field and J/^'s velocity V in the lab frame. Our results are summa- 
rized in Fig. 4. In a strong magnetic field such as the one expected to be produced at the LHC, a 
J/ip moving with P± > 9 GeV in the reaction plane is expected to dissociate because of magnetic 
field even in vacuum. At finite temperature, when binding energy decreases, dissociation becomes 
prominent at lower transverse momenta. The magnetic field has no influence on J/ip dissociation 
when it moves perpendicularly to the reaction plane (i.e. parallel to the field) because the corre- 
sponding electric field vanishes in the comoving frame. Non-trivial azimuthal angle dependence of 
w may be an important source of azimuthal anisotropy of J/ip production in heavy ion collisions. 

In Sec. V we discuss the phenomenological significance of our results. 

Before continuing with the main part of the paper, we list here some of our notational definitions 
in order to avoid confusion: V and P are velocity and momentum of quarkonium in the lab frame; 
M is its mass; p is the momentum of quark or anti-quark in the comoving frame; m is its mass; So 
is the magnetic field in the lab frame, E and B are electric and magnetic fields in the comoving 
frame; jl is the quarkonium Lorentz factor; and 7 is a parameter defined in (21). We use Gauss 
units throughout the paper; note that expressions eB, eE and eBo are the same in Gauss and 
Lorentz-Heaviside units. 

II. QUARKONIUM IONIZATION RATE 
A. Comoving frame 

Consider a quarkonium traveling with velocity V in constant magnetic field -Bo. Let B and 
E be magnetic and electric fields in the comoving frame, and let subscripts || and _L denote field 
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components parallel and perpendicular to V correspondingly. Then, 

£||=0, E ± = lL VxB , (3a) 
S|| = — B ± = 7L — 2 , (3b) 

where 7l = (1 — V 2 )~ 1 / 2 . Clearly, in the comoving frame B ■ E = 0. If quarkonium travels at 
angle <p with respect to the magnetic field in the laboratory frame, then 



5 = 5o^cos 2 <Ml- 7 2) + 7 2, E = B ~f L V sin <j). (4) 

We choose z and y axes of the comoving frame such that B = Bz and E = Ey. A convenient gauge 
choice is A = —Byx and (p = —Ey. For a future reference we also define a useful dimensionless 
parameter p [43] 

E 7 L ysin0 

5 ^COs2 0(l- 7 2) +7 2 

Note, that (i) < p < 1 because B 2 — E 2 = B^ > and (ii) when quarkonium moves perpendicu- 
larly to the magnetic field Bq, p = V. 



B. WKB method 



It is natural to study quarkonium ionization in the comoving frame. As explained in the 
Introduction, ionization is quantum tunneling through the potential barrier caused by the electric 
field E. In this subsection we employ the quasi-classical, a.k.a. WKB, approximation to calculate 
the quarkonium decay probability w. For the gauge choice specified in Sec. II A quark energy eo 
(eo < m) in electromagnetic field can be written as 

eo = \Jm 2 + (p- eA) 2 + eip = ^Jm 2 + (p x + eBy) 2 + p 2 + p 2 - eEy . (6) 

In terms of eo, quarkonium binding energy is e& = m — eo. To simplify notations, we will set p z = 0, 
because the quark moves constant momentum along the direction of magnetic field. 

The effective potential U(y) = eo(y) — \Jm 2 +p 2 corresponding to (6) is plotted in Fig. 1. Wc 
can see that the tunneling probability is finite only if E > 0. It is largest when B = 0. It has 
been already noted before in [41, 43, 44] that the effect of the magnetic field is to stabilize the 
bound state. In spite of the linearly rising potential (at B > E) tunneling probability is finite as 
the result of rearrangement of the QED vacuum in electric field. 
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FIG. 1: Effective potential U(y) — ^Jm 2 + (p x + eBy) 2 + P 2 ,— eEy — ^Jm 2 + p 2 + p 2 for p y — 7 p x = 771/6, 
B = m (except the blue line where B = 0). The width of the potential barrier decreases with E and 
increases with B. 1 — eo corresponds to the binding energy in units of m. 

Ionization probability of quarkonium equals its tunneling probability through the potential 
barrier. The later is given by the transmission coefficient 

-2/o' 1 V /Z ^^: 



w 



-f 



(7) 



In the non-relativistic approximation one can also calculate the pre-exponential factor, which ap- 
pears due to the deviation of the quark wave function from the quasi-classical approximation. Such 
a calculation requires matching quark wave functions inside and outside the barrier [45] . To deter- 
mine the pre-exponential factor in the relativistic case one needs to solve the relativistic two-body 
problem, which is analytically challenging [46]. Fortunately, as we argue later, one does not need 
to know the pre-factor to make reliable estimates of the quarkonium dissociation energy. 

We now proceed with the calculation of function /. Since B > E Eq. (6) can be written as 



Pu 



-e\B 2 - E 2 )(y - yi )(y - y 2 ) , 



where 



2/1,2 



e E - p x B T V(e E - Px B) 2 - (ff 2 - E^)(-e 2 + m 2 +g) 

e{B 2 - E 2 ) 

Define dimensionless variables eo = £o/m and q = p x /m- Integration in (7) gives: 
/ _ V-4 + 1 + Q 2 (^oE - qB) 



(8) 



(9) 



e(B 2 - E 2 ) 

(e E - qB) 2 - (B 2 - E 2 )(-e 2 + 1 + q 2 ) f e E - qB + yW~ E 2 )(-e 2 + 1 + f) 
e(B 2 - E 2 f/ 2 11 \ ^(epE - qB) 2 - (B 2 - E 2 )(e 2 + 1 + q 2 ) 



(10) 
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For different q's w = gives the corresponding ionization probabilities. The largest probability 
corresponds to smallest /, which occurs at momentum q m determined by equation [41] 

^M = . (ii) 

dq m 

Using (10) and parameter p defined in (5) we find 

pjep - pq m ) ln / epp - Qm + \/j - P 2 V~4 + 1 + qL \ = y/~~4 + 1 + 9m (u) 
!"P 2 \ " V(to-pqm) 2 -l + P 2 J 
This is an implicit equation for the extremal momentum q m = q m (^o,p)- Substituting q m into (10) 
one obtains / = f(eo,p), which by means of (7) yields the ionization probability. The quasi- classical 
approximation that we employed in this section is valid inasmuch as f(q m ) 3> 1. 

In order to compare with the results obtained in [41] using the imaginary time method, we can 
re-write Eq. (12) in terms of an auxiliary parameter To as 

T^^^ + 't^, ( 13a) 

tanhro e - pq m Mou , 
=p . (13b) 

to eoP - qm 

Taking advantage of these equations, Eq. (10) can be cast into a more compact form 

fm = — ^7==[1 - e (e - q m p)] , (14) 
eEyJl - p 2 

where we denoted f m = f(q m ). This agrees with results of Ref. [41]. Numerical solution of 
(13a), (13b), (14) is shown in Fig. 2. 



C. Special case: Crossed fields 

An important limiting case is crossed fields E = B. Since also BIB, see Sec. II A, both 
field invariants vanish. Nevertheless, quarkonium ionization probability is finite [41]. This limit is 
obtained by taking p — > 1 in the equations from the previous section. Employing (13a) and (13b) 
we get the following condition for extremum 

e 2 - 1 + 2ql - 3e q m = , (15) 



with the solution 



Substituting into (14) produces 



Qm = \ (3e - • (16) 

, _2m 2 {-el + l + glf 2 
3eE e - q m 



III. NON-RELATIVISTIC APPROXIMATION 



A very useful approximation of the relativistic formulas derived in the previous section is the 
non-relativistic limit because (i) it provides a very good numerical estimate, see Fig. 2, (ii) it allows 
us to eliminate the parametric dependence in (10), (12) and write f(q m ) explicitly in terms of p 
and eo, and (hi) spin effects can be accounted for. 



A. Arbitrary binding 

Motion of a particle can be treated non-relativistically if its momentum is much less than its 
mass. In such a case £o ~ m or = m — Eq m. Additionally, motion of a charged particle in 
electromagnetic field is non-relativistic if E <C B. Indeed, the average velocity of a non-relativistic 
particle is of order v ~ E/B = p. Thus, the non-relativistic limit is obtained by taking the limits 
Q> = £b/ m ^ 1 an d p *C 1- In these limits the extremum conditions (13a), (13b) reduce to 

to = ^ + ql , (18a) 
P 

tanhro = 

to p- q m 

Out of two solution to (18a) we pick the following one 



q m = -^T V-2e fe . (19) 
The sign of q m is fixed using (18b) by noticing that tanhro/ro < 1. Eliminating q m gives: 

^-(rocothro-l)^! 2 , (20) 

where 

(2D 

7 is analogous to the adiabaticity parameter of Keldysh [47] . Taking the non-relativistic limit of 
(14) and using (19) yields 

2m\2e b )^ 

fm = — ^ — g{i) , (22) 
where g(j) is the Keldysh function [47] 
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In Fig. 2 we show the dimensionless ratio f m eE/m 2 as a function of the binding energy e b (in 
units of m) for several values of p. The vacuum binding energy of J/ip corresponds to e b = 0.68. 
We observe an excellent agreement between the full relativistic calculation and the non-relativistic 
approximation. At p = 0.9 and e b = 0.68 the difference between the two lines is 10% and can be 
further improved by considering higher order corrections to f m [44]. 




0.0 0.5 1.0 1.5 2.0 



FIG. 2: Dimensionless function f m eE/m 2 versus tb for different values of p. The solid line is the full 
relativistic calculation, the dashed line is the non-relativistic approximation. J/tp binding energy in vacuum 
corresponds to e& = 0.68. 



B. Weak binding 

Of special interest is the limit of weak binding 7 < 1, i.e. \/2e b <C p. Expanding (20) at small 
7 and ro yields 

to = 7 (l + ^7 2 ) (24) 
and substituting into (23) and subsequently into (22) yields 

fm=~(^f /2 - (25) 

Hence, the quarkonium dissociation probability reads [45] 

W = eXP {-3^^}- (26) 
Since the quasi-classical approximation employed in this paper is valid if f(q m ) 3> 1, it follows that 
the binding energy must satisfy 

< e b < p 2 m . (27) 



m 



1/3 
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Note also that we work in the approximation of the short-range binding potential meaning that 
\[2e h < l/(mR), see Sec. I. 

C. Strong binding 

In the limit 7 > 1, (20) and (23) imply that 

ro = ^, 5(7) = y- (28) 



Substituting (28) into (22) we derive 



Thus, quarkonium dissociation probability in the case of strong binding is 



This formula is valid when 



p 2 m, y^eEp < e b < 1/R (31) 

D. Contribution of quark spin 

So far we have neglected the contribution of quark spin. In order to take into account the effect 
of spin interaction with the external field, we can use squared Dirac equation for a bi-spinor ip: 

[(e - eipf - (p - eA) 2 - m 2 + eS • B - iea ■ E] tp = , (32) 

where 



a = (33) 





Operators S z and a y do not commute. Therefore, in order to apply the WKB method for calcu- 
lation of the ionization probability one actually needs to square (32), which leads to a differential 
equation of the fourth order in derivatives. The problem becomes more tractable in the non- 
relativistic case and for crossed fields. Spin effects in crossed fields were discussed in [48] . 
With quark spin taken into account, the non-relativistic version of (6) becomes: 

^- [{p x + eByf + p 2 y ] -eEy-^s-B = -e b , (34) 
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and hence 



Py = 2m ( 



e b + ^s-B + eEy) - (jp x + eByf , 
s / 



(35) 



where [i is the quark magnetic moment and s is the projection of spin in the direction of the 
magnetic field. For a point quark, \i = hb = The effect of quark spin on quarkonium 

dissociation probability can be taken into account by replacing — > e' b = — ■ B in formulas 
for f m . With this replacement, all results of this section apply to a particle with spin. Note that 
effective binding energy e' h decreases if spin is parallel to the magnetic field and increases if it is 
antiparallel. In particular, in the case of weak binding 



Since the non-relativistic limit provide a good approximation of the full relativistic formulas, we 
will implement the quark spin dependence using the non-relativistic prescription. 



Now, we would like to use the results of the preceding sections to analyze the dissociation 
probability w of quarkonium as a function of its binding energy e&. We suppose that quarkonium 
moves with velocity V at angle 4> with respect to the external magnetic field Bq. Different binding 
energies may correspond to either different quarkonium states or to the same quarkonium state at 
different temperatures. 

Dissociation probability is given by (7) as w = exp{ — / m (e&)}. At f m {^b) 3> 1, quarkonium is 
stable with respect to Lorentz ionization because w <^ 1. The dissociation become probable only 
at such values of binding energy e& = that f m {ed) ~ 1. In view of the steep dependence of / m (ej,) 
on e&, the probability w(eb) is in fact close to the step function: it turns one when e& < and is 
zero otherwise. Consequently, we introduce the quarkonium dissociation energy in magnetic field 
£d = med as the solution to the equation 



where we took into account the contribution of quark spin, see Sec. HID. The negative sign in 
front of the spin term in (37) is chosen to give the highest dissociation energy. (Which of two 
particles core has higher probability to tunnel depends on the direction of the magnetic field). 
The solution of (37) depends on electric and magnetic fields in the comoving frame, which in turn 




(36) 



IV. DISSOCIATION ENERGY OF J/V> 




(37) 
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depend on the magnetic field in the lab frame Bq and quarkinium kinematics through the Lorentz 
transformation formulas (4). In Fig. 3 we show function / m (e&) for different velocities V of J/tp, 
assuming it moves perpendicularly to the magnetic field. In vacuum, M = 3.1 GeV, m = 1.87 GeV, 
so that £b = 2m — M = 0.64 GeV. The corresponding dimensionless parameter ej> varies in the 
interval < e b < 0.68. 




FIG. 3: f m versus for different values of J/tf) velocity V. Dissociation probability is w = exp{— f m }- 
Magnetic field in the lab frame is taken to be eBo = 15mJ. J/t/j moves perpendicularly to the field (i.e. 
in the reaction plane). Solid lines correspond to the full relativistic calculation, dashed lines to the non- 
relativistic approximation. J/ip binding energy in vacuum corresponds to t\, = 0.68. 

The numerical solution to Eq. (37) is plotted in Fig. 4. We see that the dissociation energy of 
J/ip increases with its velocity and is larger in a strong magnetic field. At eBo = 15m^, which 
according to estimates may be reached in heavy-ion collisions at the LHC [30, 31], the dissociation 
energy equals the vacuum binding energy at V = 0.94. This implies that most of J/ip'a moving in 
the reaction plane with P± > 9 GeV will dissociate in magnetic field even in vacuum. 

In a thermal medium, exhibits strong dependence on temperature T. At T = 0, equals its 
vacuum value. As T increases e b drops until it vanishes at some To near T c . The precise value of 
To as well as the functional form of £b(T) is model dependent, see e.g. [18, 19]. To illustrate the 
temperature effect, we indicated in Fig. 4 half of the J/ip binding energy £&/2 and its quarter £&/4. 
Suppose that plasma temperature is such that the binding energy is half of that in vacuum. Then, 
at eBo = 15m^, all J/^'s with V > 0.71 (P± > 3.1 GeV) will dissociate, while at cBq = this 
occurs at V > 0.92 {P ± > 7.3 GeV). 

So far in this section we have discussed the case of quarkonium moving in the reaction plane, 
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FIG. 4: Dissociation energy of J/tp as a function of its velocity V in different magnetic fields ESq. We 
assumed that J/tp moves in the reaction plane. Vacuum dissociation energy is zero. The solid horizontal 
line indicates the vacuum binding energy. Dashed horizontal lines correspond to binding energy at certain 
temperatures. 

i.e. at (p = 7r/2. Since electric field in the comoving frame is proportional to sino\ see (4), the 
dissociation probability decreases at oblique angles. J/ip's moving parallel to the magnetic field 
Bq, i.e. at 4> = 0, are not affected by electromagnetic field at all. This non-trivial azimuthal angle 
dependence translates into azimuthal anisotropy of produced J/V's, which we plan to address in a 
separate paper. 

The size of the spin contribution is determined by the dimensionless ratio £; = eB /2et,m 2 , see 
(36). Its magnitude for J/tjj moving with a non-relativistic velocity perpendicularly to the magnetic 
field Bq is £ « (ei?o/240m^)(0.68/e;,). Therefore, even in the strongest field expected at the LHC, 
eBo = 15m 2 , the effect of spin is quite small at T = 0. £; increases as the temperature of plasma 
and/or velocity of J/ip increases owing to decrease of the binding energy from its vacuum value 
€b = 0.68 down to zero. When £ ~ 1 non-relativistic treatment of spin breaks down. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we chose the quarkonium proper frame to calculate the dissociation probability. 
In this frame the potential energy of, say, antiquark (with e < 0) is a sum of its potential energy in 
the binding potential and its energy in the electric field —eEy, where y is the electric field direction. 
Since \e\Ey becomes large and negative at large and negative y (far away from the bound state) 
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and because the quarkonium potential has finite radius, this region opens up for the motion of the 
antiquark. Thus there is a quantum mechanical probability to tunnel through the potential barrier 
formed on one side by the vanishing quarkonium potential and on the other by increasing absolute 
value of the antiquark energy in electric field. Of course the total energy of the antiquark (not 
counting its mass) is negative after tunneling. However, its kinetic energy grows proportionally to 
eEy as it goes away. By picking up a light quark out of vacuum it can hadronize into a D-meson. 

If we now go to the reference frame where E = and there is only magnetic field B (we can 
always do so since E < B), then the entire process looks quite different. An energetic quarkonium 
travels in external magnetic field and decays into quark-antiquark pair that can late hadronize into 
.D-mesons. This happens in spite of the fact that J/ip mass is smaller than masses of two D-mesons 
due to additional momentum eA supplied by the magnetic field. Similarly a photon can decay into 
electron-positron pair in external magnetic field. 

In the previous section, we observed that J/ip dissociation energy increases with magnetic field 
strength and quarkonium momentum. As a consequence, quarkonia dissociate at lower temperature 
than one would have expected based on calculations neglecting magnetic field. Fig. 4 indicates that 
in heavy-ion collisions at the LHC, J/^'s moving with P± > 9 GeV in the reaction plane dissociate 
with probability of order unity even in vacuum, where we rely on the magnetic field estimates of 
[30] . This effect is more pronounced at finite temperature owing to the decrease in binding energy 
and, as seen in Fig. 4, it occurs in weaker fields. 

Although magnetic fields in pp and pA collisions are much weaker than in AA collisions, they are 
still strong enough to cause J/ip dissociation at sufficiently high momenta P±. A truly spectacular 
feature of such process would be J/ip decay into two heavier D-mesons. Validity of the quasi-static 
approximation can perhaps be advocated for a pre-hadronic cc state owing to its short production 
time. This intriguing possibility will be addressed elsewhere. 

The effect of J ftp dissociation in a magnetic field vanishes in the direction parallel to the 
magnetic field, i.e. perpendicular to the reaction plane. Therefore, J/tp dissociation gives negative 
contribution to the total azimuthal asymmetry coefficient V2- The size of this contribution is a 
subject of our ongoing investigation. 

Finally, we would like to mention a possibility of J/tp dissociation directly by electric fields in- 
duced in plasma by the Chiral Magnetic Effect [30].' We plan to address this interesting possibility 
in a forthcoming publication. 



This observation is due to D. Kharzeev. 
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